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Abstract 

A leading-order analysis of E866/NuSea and NA3 Drell-Yan data in nuclei is car- 
ried out. At Fermilab energy, the large uncertainties in the amount of sea quark 
shadowing prohibit clarifying the origin of the nuclear dependence observed exper- 
imentally. On the other hand, the small shadowing contribution to the Drell-Yan 
process in tt~-A collisions at SPS allows one to set tight constraints on the en- 
ergy loss of fast quarks in nuclear matter. We find the transport coefficient to be 
q = 0.24 lb 0.18 GeV/fm^ that corresponds to a mean energy loss per unit length 
-dE/dz = 0.20 ± 0.15 GeV/fm for Eg > 50 GeV quarks in a large {A w 200) 
nucleus. 
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1 Introduction 



Energy loss of hard partons in hot QCD matter is expected to be large 0. 
Consequently, it has been suggested that the depletion of high p± jets (jet 
quenching) in heavy ion collisions may be a signal for quark-gluon plasma 
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formation 0, Recently, a lot of excitement was created as PHENIX data 
revealed that high p±_ hadron spectra were found to be substantially suppressed 
in the most central Au-Au collisions at RHIC with respect to the extrapolation 
from p-p data f^. Multiple scattering of a high energy parton traversing a 
large nucleus ("cold" QCD matter) has been studied similarly by Baier et al. 
(BDMPS) in Ref. in which a numerical estimate for the expected parton 
mean energy loss per unit length —dE/dz is given. 

The Drell-Yan mechanism is a process particularly suited for the study of 
quark energy loss in nuclei as the lepton pair does not strongly interact with 
the surrounding medium. Furthermore, new p-A data recently became avail- 
able from the E866/NuSea experiment at Fermilab 0]. Subsequently, two re- 
cent attempts to extract the quantity —dE/dz from these data have been car- 
ried out P, 1^ . However, their results do not agree as the amount of sea quark 
shadowing assumed in both analysis strongly differs. As we shall see later, 
the poorly known shadowing corrections at Fermilab energy indeed makes a 
model-independent extraction of quark energy loss unlikely. 

In this letter, we discuss constraints on —dE/dz from the analysis of both 
E866/NuSea [^] and NA3 [^] Drell-Yan dimuon data in hadron nucleus reac- 
tions. The procedure followed is detailed in Section 2 after having given the 
leading-order Drell-Yan production cross section in nuclei. The results given in 
Section 3 are discussed and compared with previous studies in the last section. 



2 Nuclear dependence of Drell-Yan production 



Leading order production cross section 

To leading order (LO) in perturbation theory, the Drell-Yan (DY) process de- 
scribes dilepton production from quark-antiquark annihilation. The invariant 
mass M of the lepton pair is set by the center-of-mass energy of the qq colli- 
sion = (xiX2s)^/^ where Xi (resp. X2) is the momentum fraction carried 
by the beam (resp. target) parton and ^/s is the center-of-mass energy of the 
hadronic collision. The differential partonic cross section qq — >■ has been 
computed to leading order in § and is given by 



da Stto"^ 2 



dM 9M 



eiS{s-M'). (1) 



The hadronic DY production cross section is then obtained from the convo- 
lution of the partonic cross section (1) with the quark distributions in the 
beam and in the target hadron, evaluated at Xi and X2 respectively, and at a 
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factorization scale /i^ = M^, i.e. 
da{hh') Svra^ 



dxidM 9M xis g 



E 4 (/.'(^O/f (^2) + (x,)) (2) 



where the sum is carried out over the hght quark sector q = u,d,s and 
X2 = M"^ /xis after integration over the delta function[^. Looking at the Xi 
dependence of DY production, we shall rather deal in the following with the 
single differential cross section 

da{hh') 8na^ ^ ^ r dM 



dxi 9xi s 



T.4j^ (/^(^l)/,-'(^2) + fi{x,)f^'{x2)) (3) 



where the integration over the dilepton mass is performed in the range between 
the cc and the bb resonances. The LO cross section (3) proves to describe 
the whole trend of DY data within a so-called ~ 2 factor which might 
be attributed to large next-to-leading order (NLO) corrections. Since we are 
primarily interested here in the nuclear dependence of Drell-Yan production, 
one may reasonably expect these higher-order corrections to cancel in the 
production ratio 

h,., . B (daihA)\ (daihB)\~^ 



in a heavy (A) over a light (B) nucleus. We shall therefore restrict ourselves 
to a LO analysis throughout this paper. In the absence of nuclear effects, the 
production cross section in hadron-nucleus reactions a{hA) appearing in (4) 
is given by 



^(/,'(^i)/|(^2) + /^(xi)/,^(a;2)) (5) 
+ {A - Z) {fiix,)f^{x2) + 

after separating into terms involving protons and neutrons in the target nu- 
cleus. Isospin effects will remain small as long as the parton densities in the 
proton and in the neutron do not differ strongly ff{x2,f^^) ~ /"(a;2,/i^), that 
is, at small X2 = M'^/xis <^ 1. In this particular case, DY production scales 




^ In the following, parton distributions will always be evaluated at the hard scale 
/i^ = M^. For simplicity, we shall drop the explicit dependence and use fi{x) = 
fi{x^lj? = M^) in the notations. 
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with the atomic mass number A, hence the nuclear production ratio (4) is 
equal to one. 

Beyond isospin corrections, several nuclear effects — such as shadowing or 
parton energy loss — might affect the Drell-Yan process and lead to an un- 
usual A-dependence in dilepton production. Let us discuss now how these 
mechanisms modify the nuclear production cross section (5). 



Shadowing 

Nuclear deep inelastic scattering (DIS) data indicate that parton distributions 
in nuclei differ significantly from those in a proton [|Ty]. In particular, a sig- 
nificant depletion ("shadowing") of high momentum (0.3 < X2 < 0.7) as well 
as low momentum (x2 < 0.05) partons in large nuclei has been reported ex- 
perimentally. The origin of nuclear shadowing is still rather unclear. It may 
be attributed to the multiple scattering of the struck quark with the target, 
which can be removed into effective nuclear parton distribution functions ac- 
cording to the QCD factorization theorem |TT|. Therefore, Drell-Yan dilepton 
production off nuclei can be estimated using the nuclear densities ff^^ (resp. 
f^^^) instead of the "free" parton distributions /f (resp. /f) in Eq. (5). In 
the following numerical applications, we will assume that the nuclear parton 
distributions /f^"^ factorize into a nuclear contribution Rf and the parton 
distribution in a proton /f , i.e., 

f!/^{x,f^') = Rt{x,fj.')xf[{x,ix') (6) 



where the function Rf{x, /x^) has been parameterized by Eskola, Kolhinen, and 
Salgado (EKS98) from a leading-order DGLAP analysis of DIS data 0. It 
is worth pointing out that E772 Drell-Yan data have also been taken into 
account in their fitting procedure to further constrain sea quark shadowing 
in the intermediate X2 range {x2 ~ 0.1) [|12|. We shall come back to this 
observation when discussing results in Sections 3 and 4. 



Parton energy loss 

The quark (antiquark) from the projectile may scatter through the nucleus 
before the hard qq annihilation process occurs. The medium induced gluon 
emission from the incoming parton with energy Ep leads to a radiative parton 
energy loss e. In the BDMPS approach, the distribution D{e) in the energy 
loss is characterized by a typical energy scale uJc proportional to the square of 
the length L of traversed nuclear matter 

u^c = lqL'. (7) 
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The so-called "transport coefficient" q relates the p±_ broadening of the par- 
ton to the length L 0]. It is expected to depend on both the small x gluon 
distribution and the density p of scattering centers of the medium (here, 
p = 0.15 fm~'^). Assuming the hard process to take place uniformly in the 
nucleus, the length L is proportional to the nuclear radius, L = 3/4 i?. 



Neglecting interference effects in the multiple gluon radiation, Baier et al. give 
the distribution D{e) a simple integral representation For the applications 
to come, the distribution D{e) has been computed numerically in the soft gluon 
approximation, e ~ c^c Ep, from Eq. (18) of Ref. [Q. The mean BDMPS 
energy loss AE of the incoming quark is given by M 



-AE = JdeeD{e) = ^asCRUJcOcL^ (8) 

with Cr being the color charge of the parton {Cr = 4/3 for quarks) and 
as = 1/2 the strong coupling constant. In the following, we shall therefore 
write the mean energy loss per unit length as 



where 5 is a free parameter simply related to the transport coefficient q through 
(7) and (8). 

This multiple scattering effect shifts the quark (antiquark) momentum fraction 
from Xi + Axi(e) to Xi at the point of fusion, with 

Ax.ie) = ^ (10) 



and where is the projectile hadron energy in the nucleus rest frame[^. 
Consequently, the parton densities fl^{x) have to be evaluated at (xi -|- Axi) 
in the nuclear production cross section (5). Because of the steep behavior of 
the valence quark distributions at large Xi (e.g., ~ (1 — Xi)^~^ in a proton), 
even a small shift Axi may substantially suppress Drell-Yan production in a 
large nucleus as compared to a hght one. We further note that the larger xi, 
the stronger the suppression Uy{xi + Azi)/m^(xi). 



Analysis of Drell-Yan data 

In the previous sections we have stressed that nuclear mechanisms affect the 
Drell-Yan process. In the most general case, the production cross section in 

^ Let us drop in the following the explicit e dependence of Axi for clarity. 
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hadron-nucleus reactions will read 



MkA)_S..'^^,jiMr^^^^^ (11) 



dxi 9xis ^"'^ J M 
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+ /\xr)fl'\x^) + {A- Z)f^^{x, + Axi)/;/^(x2) 

taking both nuclear shadowing and parton energy loss effects into account. 
At large xi, the restricted phase space (e < (1 — xi) Eh) makes the effects of 
quenching even more pronounced. In the limit of no shadowing {ff^"^ = ff) 
and vanishing energy loss (Axi = 0), one retrieves the usual production cross 
section (5). 

The aim of our study is to investigate whether the above mentioned effects 
manifest themselves in available data, and, if so, to possibly disentangle shad- 
owing from energy loss contributions to the nuclear dependence of Drell-Yan 
production. In particular, it would be most interesting to set some constraints 
on the amount of parton energy loss in nuclear matter. To achieve such a 
goal, a close comparison between data and theory has been carried out. Using 
Eq. (11), Drell-Yan production is computed with the following four options: 

(i) /f = /f , Axi = 

(ii) /f ^ /f , Axi = 

(iii) /f = /f , Axi ^ 

(iv) ff^ ^ /f , Axi ^ 0. 

No nuclear effect is assumed in the first set (i), while shadowing and energy 
loss corrections are considered in turn ((ii) and (iii) respectively). Both effects 
are then combined in the last case (iv). Shadowing corrections were taken from 
the EKS98 parameterization (Eq. (6)) and the momentum fraction shift Axi is 
given by Eq. (10) with 6 (and hence the coefficient q) kept as a free parameter 
fitted to the data. We made use of the MRST LO parton distributions in 
a proton /f |jl6[ and their similar study in the pion f[ |T^. The parton 



distributions in the neutron /" as well as shadowing corrections fl"^"^/ ff' are 
given by the proton (nuclear) distributions with the usual assumptions: = 

rf", dP = u", uP = J", dP = u", and sP = s"". 

The theoretical calculations were then confronted separately with two sets of 
data. First, the E866/NuSea collaboration reported recently on high-statistics 
measurements of Drell-Yan dimuon production in proton-nucleus (Be, Fe, W) 
collisions using the 800 GeV proton beam at Fermilab P . They extracted both 
production ratios -R(Fe/Be, xi) and i?(W/Be, xi) over a large kinematic accep- 
tance (0.28 < xi < 0.84) and on the 4.0 < M < 8.4 GeV mass range. The sec- 
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ond data set corresponds to older measurements from the NA3 collaboration 
of DY production in pion induced reactions at beam energies i?^- = 150 GeV 
and Et,- = 280 GeV 0. Although the statistics are somehow limited, the ratio 
of production off hydrogen over platinum targets, (p/Pt), was extracted 
at both energies and up to xi ~ 0.9 in the mass range 4.1 < M < 8.5 GeV. 



3 Results 

E866/NuSea data 

Let us first present the results from the fit to the Fermilab data which consist of 
both ratios i?(Fe/Be, xi) and i?(W/Be, xi) (7 data points each) in p(800 GeV)- 
A collisions. The agreement between data and theory is summarized in Table 1 
where the per number of degrees of freedom is given. 

First, it is clear from Table 1 that the E866/NuSea data exhibit significant 
nuclear effects from the large /ndf = 4.49 when neither shadowing nor en- 
ergy loss is taken into account. Even though nuclear shadowing and/or energy 
loss gives an excellent description {x^/ndf ^ 0.5), these data do not allow one 
to pin down one or the other case from the constant x^/nd/ in each separate 
scenario. Perhaps more interesting is the amount of energy loss required to 
describe E866/NuSea measurements. When EKS98 shadowing is not included 
(iii), it turns out that a large 6 = 3.5 GeV/fm is required which corresponds, 
using Eq. (9), to an energy loss per unit length of dE/dz = 1.75 GeV/fm in 
a large nucleus (L ^ 5 fm). This result turns out to be close to (although 
well smaller than) the recently fixed dE/dz ^2.7 GeV/fm by Johnson et 
al. from the Xi and M dependence of E772 and E866/NuSea data. On the 
contrary, no significant energy loss (5 = 0.1 GeV/fm) is found when EKS98 
shadowing is included in the calculations (iv), which confirms previous results 
of Ref. P]. However, one should keep in mind that E772 Drell-Yan measure- 
ments, taken in the same kinematic range as E866/NuSea, have been used 
to constrain the EKS98 parameterization. This may therefore explain why 
E866/NuSea data (consistent with E722 results) are well reproduced assum- 
ing EKS98 shadowing only {x^/ndf = 0.51, (ii)). This lack of consistency 
thus weakens our confidence in a vanishing energy loss in Fermilab Drell-Yan 
data. Therefore, quark energy loss appears to strongly depend on the initial 
assumptions, as illustrated in Figure 1 (left) where the /ndf is plotted as a 
function of 5. Moreover, Figure 1 shows that the two fitted values are rather 
well constrained from the deep minima. To quantify the error on the parame- 
ter 5, we display between brackets in Table 1 the upper limits 5 + where 

^ The number of degrees of freedom is 14 and 15 depending on whether the energy 
loss coefficient 5 is taken or not as a free parameter. 
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the one standard- deviation error AS is given by the deviation of by one 
unit from its minimum. 

Consequently, E866/NuSea data do not permit conclusions on the very origin 
of the observed nuclear dependence because of the present uncertainties in the 
amount of sea quark shadowing. This prevents us from setting constraints on 
quark energy loss in nuclear matter. Let us now turn to the NA3 DY dimuon 
production in pion nucleus reactions. 





(i) 


(ii) 


(iii) 


(iv) 


E866/NuSea 


6 (GeV/fm) 






3.5 (4.0) 


0.1 (0.6) 


X^/ndf 


4.49 


0.51 


0.52 


0.54 


NA3 


5 (GeV/fm) 






0.3 (0.7) 


0.5 (1.0) 


X^/ndf 


0.38 


0.43 


0.39 


0.38 


both 


6 (GeV/fm) 






1.6 (1.8) 


0.3 (0.7) 


X^/ndf 


2.32 


0.48 


1.34 


0.47 



Table 1 



Results from a fit to E866/NuSea to NA3 data sets for the various ansatz assumed in 
the calculations (see text). The 1 a upper limits 5 + /S.5 are given between brackets. 



NA3 data 

Unlike the DY data taken at Fermilab, Drell-Yan production in vr" induced 
reactions at E^,- = 150 and 280 GeV should not be spoiled by large shadowing 
corrections. The reason for this is twofold. First, the mean momentum fraction 
X2 probed in the NA3 measurements turns out to be much larger (0.06 < 
(X2) < 0.3) than at Fermilab (0.02 < (^2) < 0.06) because of the smaller 
incident energy. In this intermediate X2 range (x2 ~ 0.1), DY should only be 
slightly affected by quark (anti)shadowing |[T^. Furthermore, the Drell-Yan 
process in these reactions is dominated by the annihilation of valence quarks 
for which shadowing is well constrained from DIS measurements only. 

The resuhs from the fit to the ratio R""' {p/Pt, xi) at E^- = 150 and 280 GeV 
(resp. 8 and 9 data points) are displayed in Table 1. First, Table 1 shows 
that the NA3 data are well accounted for without invoking any nuclear effect 
{x'^/ndf = 0.38), as already pointed out in The agreement can be seen 
in Figure 2 (left) where the measurements at 150 GeV are compared to the 
calculations ((i), so/z(i)0. Furthermore, we note that shadowing effects only 

^ Notice that data are consistent with the expected (p/Pt) = 2A/{Z + A) ^ 
1.43 assuming only valence- valence fusion process in Eq. (5) with u'^{x2) « 2uy{x2) 
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marginally affect DY production ((ii), dashed) as expected. As a consequence, 
we anticipate that energy loss can be fixed from tliese data witliout ambiguity 
anymore. 

^ 5 I 

■a 
a 

4-5 ^ 

4 : 
3.5 : 

3 : \ 
2.5 : 

2 : 

1.5 : ""X' 

1 : \,.. 
0.5 

"o 0.5 1 1.5 2 2.5 3 3.5 4 "o 0.5 1 1.5 2 2.5 3 3.5 4 

5 (GeV/fm) 8 (GeV/fm) 

Fig. 1. jndf between E866/NuSea (/e/t) and NA3 (right) data and tfieoretical 
calculations without ((hi), dotted) and with ((iv), dash-dotted) shadowing corrections 
as a function of the energy loss coefficient 6. The jndj for vanishing energy 
loss ((5 = 0) is also shown without ((i), solid) and with ((ii), dashed) shadowing 
contribution. 

Within both ansatz (iii) and (iv), NA3 data reveal that quark energy loss is 
small, from 5 = 0.3 GeV/fm up to 0.5 GeV/fm when shadowing is taken into 
account. What is more, the small one standard-deviation upper limits (0.7 
and 1.0 GeV/fm) clearly indicate that DY measurements at SPS — despite 
the large error bars — put stringent constraints on the maximal quark energy 
loss in nuclear matter. The strong disagreement between NA3 data and theory 
for too large energy loss coefficients S is shown in Figure 1 (right) where the 
X^/ndf is plotted. In particular, they allow to exclude the huge energy loss 
6 — 3.5 ± 0.5 GeV/fm extracted from the E866/NuSea DY data assuming 
no shadowing effects (cf. Table 1, (iii)). To get a feeling for the origin of 
such tight constraints, we plot in Figure 2 (right) the theoretical predictions 
for 5 = (solid), 1.5 (dashed) and 3 (dotted) GeV/fm in comparison to the 
150 GeV data. There, we see that the effects of quark energy loss become 
significant at large xi, leading to a fast increase of the ratio (p/Pt, xi) in 
contradistinction to the trend of the data. Even though the increase of /ndf 
comes mainly from the region of large Xi measurements, it is worth noting 
that data above Xi > 0.6 cannot be accommodated with, say, 6 = 3 GeV/fm. 
From the combined fit to NA3 data with (iv) and without (iii) shadowing, we 
found the energy loss coefficient to be 5 = 0.4 ± 0.3 GeV/fm for fast quarks 
in nuclei. We shall discuss this result in the next section. 

in this X2 range. 



■o 

-^4.5: 
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Finally, a global analysis of the E866/NuSea and NA3 overall data has been 
performed. As expected, Table 1 not only indicates that data show a significant 
nuclear dependence {x^/ndf = 2.32) but clearly demonstrate that the energy 
loss mechanism alone cannot account for it {x^/ndf = 1.34). Hence, this is 
a clue that shadowing is at the origin of the nuclear dependence observed at 
Fermilab. 



4 Discussion 



The result of our analysis is that NA3 pion-nucleus data proved more effective 
than the Fermilab precise measurements to constrain the amount of quark 
energy loss in matter which turns out to be small: 6 = 0.4 ± 0.3 GeV/fm. 
Before discussing this result and comparing it with previous studies, we would 
like to comment on limitations and uncertainties in this approach. 

First, the tight constraints on maximal quark energy loss arise mainly from 
the large xi measurements at 150 GeV (see Figure 2) as already discussed. 
To investigate the sensitivity of these data points on our final result, a simi- 
lar procedure has been performed removing, respectively, the largest (xi = 
0.95) and the three largest valued (xi > 0.7) Xi data points, leading to 
6 = 0.5 ± 0.4 GeV/fm and 6 = 0.7 ± 0.5 GeV/fm. Surprisingly, quark en- 
ergy loss remains well constrained, although a bit larger, from the whole NA3 
measurements. Furthermore, we have checked that the theoretical calculations 
depend only marginally on a specific choice for the proton parton distributions 
in Eq. (11). Indeed, similar results were found using either the leading order 
GRV LO ig and CTEQ5L |T9[ or the next-to-leading order MRST |g parton 



densities. Turning to the pion sector, it is regrettable that quark distributions 
are much less constrained from the limited available data. Let ~ (1 — x)^ 
be the valence quark distribution in the pion at large x. In Ref. |T^, the 
analysis of NAlO and E615 data lead respectively to ?7 = 1.08 ± 0.02 and 
7] = 1.15 ± 0.02. Nevertheless, going from t] = 1.06 to ?7 = 1.17 only affected 
our final result (assuming = 1.11) by a few percent at most. On the theoret- 
ical side, the Drell-Yan mechanism has been computed in the QCD-improved 
parton model to leading order in the coupling, assuming the NLO corrections 
to vanish in the production ratio R. This ansatz is justified as long as nuclear 
effects involved in both LO and NLO processes remain identical, which may 
not be true in general. Consider for instance the Compton scattering qg — >■ qj* 
process. The larger gluon energy loss {Cr = 3 in Eq. (8)) would lead to an even 
smaller quark energy loss estimate^. However, a complete calculation includ- 



^ On top of a different gluon energy loss, antishadowing may be more pronounced 
for gluons than for valence quarks in the X2 ~ 0.1 window whereas small X2 shad- 



owing in the gluon and the sea quark channel should be quite similar |12|| . 
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ing all NLO processes would certainly be needed. Concerning the computation 
of quark energy loss, the distribution D{e) has been calculated in the soft gluon 
limit Uc <^ Eg. It has been a posteriori checked that, with 6 = 0.4 GeV/fm, 
this relation is fulfilled for quark energy down to Eg 50 GeV (i.e., the small- 
est quark energy in the NA3 data). Furthermore, the Glauber approximation 
(i.e., that multiple successive quark-nucleon scatterings are independent) on 
which the BDMPS framework relies is only relevant for highly energetic quarks 
and should break down at very large energies when shadowing effects become 
large Therefore, we can assume that the use of a Glauber based approach 
to describe moderate and large Xi data at SPS energy (not too small X2) is 
sufficiently accurate. As for shadowing corrections, we believe the factorized 
form (6) assumed in Ref. |T^ to be meaningful when X2 is not too small, 
where DY production has been measured {x2 > 10~^). What is more, the use 
of the EKS98 parameterization which is well constrained from both DIS and 
DY data — although not fitted to the above NA3 measurements — appears to 
be fully justified. 



3.5 



2.5 

2 

1.5 
1 
0.5 




— no shadowing 
with shadowing 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X, 



3.5 



Q 3 
05 



2.5 
2 

1.5 
1 

0.5 




5 = GeV/ fm 
5 = 1.5 GeV/ fm 
8 = 3 GeV/ fm 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X, 



Fig. 2. NA3 ratio (p/Pt) of Drell-Yan dimuon production versus xi in 7r~(150 
GeV)-A collisions. Calculations assuming shadowing (left) and energy loss {right) 
effects are compared to the data. 

We pointed out that several attempts to extract quark energy loss in nuclear 
matter from the E772 and E866/NuSea Drell-Yan data have recently been 
made 0, 0. In Ref. 0, large energy loss effects have been claimed to be ruled 
out using E866/NuSea data corrected for shadowing and assuming various 
models for parton energy loss. Using the BDMPS approach, the one standard 
upper limit was found to be 5 + A6 = 0.46 GeV/fm, comparable to what 
is quoted here. Nevertheless, as mentioned repeatedly, their result should be 
taken with care since the shadowing parameterization is partially fitted to 
E772 data. To avoid such an inconsistency, a new analysis of both E722 and 
E866/NuSea results came out later in which nuclear shadowing was estimated 
theoretically M. Because of the small shadowing effects in their calculation, a 



11 



huge quark energy loss dE/dz = 2.7±0.4±0.5 GeV/fm proved necessary to 
describe the overall Fermilab data set, i.e., rather close to (although somewhat 
larger than) our estimate dE/dz = 1.75 ± 0.25 GeV/fm (with L = 5 fm) from 
the E866/NuSea analysis assuming no shadowing (Table 1, (iii)). Such a large 
number cannot be understood as coming from radiative energy loss only 
It has been attributed in Ref. to the interplay between the effects of gluon 
radiation together with the energy loss due to the string tension. The result 
quoted in Ref. is likely to be ruled out from a comparison with the NA3 
measurements. In particular, our present work does not show any evidence 
for quark energy loss coming from the tension of the string stretched from 
the beam parton to the nucleus. Finally, we would like to stress that the 
distribution in the induced energy loss D{e) has been employed in the present 
study whereas the previous analysis f^, |^ modeled the quenching by shifting 
the projectile quark energy by the mean energy loss AE Q. As recently pointed 
out by Baier et al., this standard modeling of the suppression is inadequate 



when the cross section sharply falls down with xi [|I^. In particular, it is 
argued that the typical quark energy loss, i.e., the loss that really contributes 
to the quenching, proves much smaller than the mean energy loss AE. 

It is interesting to note that —dE/dz = (0.4±0.3 GeV/fm) (L/lOfm) revealed 
from the NA3 data proves in excellent agreement with the BDMPS predic- 
tionQ -dE/dz = (0.4 GeV/fm) (L/lOfm) in Ref. @]. However, this apparent 
agreement is coincidental because of both the quoted error bars and the rough 
estimate in 0]. In particular, let us stress that the NA3 data are compatible 
with zero energy loss (Table 1). Nevertheless, the "agreement" is a hint that 
the origin of the fitted energy loss is radiative. This allows us to extract the 
"transport coefficient" for cold QCD matter q = 0.24 ib 0.18 GeV/fm^, which 
corresponds to a.p± broadening dp\/dz = 0.021 ±0.016 GeV^/fm in agreement 
with E772 results [|l^]. Unfortunately, such a quantity has not been measured 
by the NA3 collaboration to check the relation between radiative energy loss 
and pj_ broadening. 

To summarize, a LO analysis of Drell-Yan data in p-A and tt~-A reactions 
at Fermilab and SPS energies has been carried out. The aim was to set tight 
constraints on quark energy loss in nuclear matter in a (as much as possible) 
model-independent way. For this, multiple fits to the data have been per- 
formed under various assumptions as for the nuclear effects. At Fermilab en- 
ergy, Drell-Yan measurements probe a small X2 range where the amount of sea 
quark shadowing is only poorly known from DIS data. Consequently, we were 
unable to explore the origin of the nuclear dependence seen in the data and 
thus to extract unambiguously quark energy loss. On the other hand, nuclear 



^ This would correspond to take D{e) = (5(e — AE) in Eq. (11) 

^ The prediction in Ref. [H] has to be increased by a factor of two from elastic 



corrections originally not taken into account [15|. 
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shadowing of valence quarks in the intermediate X2 range is rather small and 
well constrained from DIS data. Therefore, DY production in n-A collisions at 
SPS energies is only sensitive to energy loss of the incoming valence antiquark. 
Prom NA3 measurements, we give the estimate —dE/dz — 0.20±0.15 GeV/fm 
for the quark mean energy loss in a Platinum nucleus, in good agreement with 
the expectation from the BDMPS perturbative approach. While smaller error 
bars would be needed to fix it more precisely, these data already allow one 
to rule out a quark mean energy loss much greater than 0.5 GeV/fm. This 
small radiative energy loss thus gives a hint that most of the Drell-Yan nu- 
clear dependence at Permilab actually comes from large shadowing corrections 
(as first assumed in the EKS98 parameterization) which may be even stronger 
at RHIC and LHC energies. 
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